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Abstract 

Results of 2005-2007 campaign of measurement of the optical turbulence 
vertical distribution above Mt. Maidanak are presented. Measurements are 
performed with the MASS (Multi-Aperture Scintillation Sensor) device which 
is widely used in similar studies during last years at several observatories across 
the world. The data analysis shows that median seeing in free atmosphere (at 
altitudes above 0.5km) is &!46 and median isoplanatic angle is 2.'47. Given a 
rather long atmospheric coherence time (about 7 ms when the seeing is good) 
such conditions are favorable for adaptive optics and interferometry in the 
visible and near-IR. 
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1 INTRODUCTION 



Studies of atmospheric optical turbulence have progressed significantly in the last decade. This 
was stimulated by new programs of site testing for extra-large telescopes. It was also realized 
that the potential of existing observatories may be further enhanced by the use of adaptive 
optics (AO) systems. 

It is well known that knowledge of only the integral turbulence on the line of sight is not 
enough for correct prediction of AO efficiency at a given site. Information on the vertical distri- 
bution of turbulence is needed for the development of a particular kind of AO (see for example 
Vernin et all Il99ll : lEoddierl . Il999l : IWilson et all . 1200.4 iFuensalida et all 120041 : ITokovinin etaL 



2003al ). It is also well known that atmospheric limits on photometric and astrometric precision 



are directly linked to th e vertical profile of turbulence ( Dravins et al. . 1997 : Kenyon et al. . 20061 : 



Shao fc Colavhal . ll992h . It is worth noting that the discussion of this issue in the 1970-80th 



was mainly related to the selection of optimal height of telescope towers. 

Astroclimate parameters of Maidanak observatory were studied for m any years. One such 
camp aign w as conducted in 1998-99, as reported by Ilyasov et al. ( 19991 ) and Ehgamberdiev 
et al (|2000h. First est i mates of the contribution of free atmosphere to seeing were made by 



(|Kornilov Tokovininl . l200lh aased on stellar scintillation analysis. It turned out that this 
contribution was quite significant, about 30%. 

The developm ent of an effective technique to study vertical turbulence distribution by 
stellar scintillations ( Tokovinin Kornilov . 2002 : Kornilov et al. . 20031 : ITokovinin et all l2003bl ) 
and of respective instrumentation allowed us to initiate in 2005 the project of turbulence moni- 
toring at Mt. Maidanak during 2-3 years, completed in 2007. Here we report its main results. 



2 MEASUREMENT METHOD AND MULTI-APERTURE 
SCINTILLATION SENSOR (MASS) 

The measurement of the altitude distribution optical turbulence is based on the fact that am- 
plitude of stellar scintillation produced by a turbulent layer depends on its height in a different 
way when observed throu gh different apert u res. 



As well described in 



Tatarskivl (|1967l ) ; iRoddierl (|l98ll ) , stellar scintillation is produced by 



phase fluctuations of a light wave passing through optically turbulent layers of atmosphere, and 
subsequent propagation. The strength and spatial spectrum of amplitude fluctuations produced 
by the atmosphere depend on the propagation length. 

Scintillation is quantified by t he sci ntillation index s 2 representing the relative variance 
of flux in some aperture. Tokovinin ( 19981 ) generalizes this concept to the case of simultaneous 
measurement of fluxes in two different apertures. The two kinds of scintillation indices - so called 
"normal" and differential - can be directly measured. On the other hand, these indices equal 
integrals of the turbulence distribution along the line of sight with some weighting functions. The 
weights are computed from the known aperture sizes and spectral energy distribution of the light. 
It should be noted that the validity of the small-perturbations method and of the Kolmogorov 
turbulence spectrum is essential here. In typical conditions of astronomical observations, both 
these assumptions are usually valid. 

Further dev elopme nt of this method is described in papers by Kornilov et al. (|2003h and 
Tokovinin et al. d2003bl ). where specific algorithms of restoration of turbulence along the line 
of sig ht are presented and various i nstrumental fac t ors a re accounted for. As dedicated studies 



show (ITokovinin &: Kornilovl . 120071 : iKornilov et all 120071 ) , the precision of turbulence strength 
representation by a set of six altitude layers is de finitely better tha n 10% . An independent 
cross-check of these results at some observatories ( Fuensalida et all 120071 : ISadibekova et al. 
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2006) was made with a SCIDAR device (Scintillation detection and ranging. lAzouit &z Vernin 
1980; iFuchs et all Il998h . 



Table 1: MASS entrance aperture diameters ( d\ — inner diameter, do — outer diameter). 



Channel 
MASS 


Segmentator mirrors 
di, mm do, mm 


Effective aperture 
di, cm do, cm 


A 




1.03 




1.60 


B 


1.05 


1.90 


1.62 


2.93 


C 


1.95 


3.60 


3.00 


5.54 


D 


3.65 


6.75 


5.62 


10.4 



The MASS device is simply a fast 4-channel photometer which measures fluctuations of 
the light flux from some reasonably bright star in four concentric ape rture s from 2 to 10 cm 
diameters. Detailed description of the device is given byKornilovet al (|2003l ). Any 5- to 6-inch 
telescope of about 3m focal distance and without central obscuration may be used as a feeding 
optics. 

A natural drawback of the MASS method is its insensitivity to turbulence near the ground. 
In order to include it for accessing complete turbulence integral along the line of sight (giving 
the seeing measu re), one has to use some other method, for example differential image motion 
monitor (DIMM. ISarazin fc E,oddierl . [l99(ih . It is worth recalling here that image quality ( seeing ) 
is defined as a full-width at half-maximum of a long-exposure stellar image created by a large 
aberration-free telescope. Seeing produced by the whole atmosphere without a contribution from 
the ground layer is referred to as free atmosphere seeing (3 f rP ^- A reliable e stima t e of this quantity 
is norm ally obtained for turbulence above h 
2003bl ). 



0.5km (jKornilov et all 120031 : iTokovinin et al 



The altitude distribution of turbulence above Mt. Maidanak was studied with help of 
so-called original MASS device of first generation which was designed and built in 2001-2002 
in the Sternberg Astronomical Institute in collaboration with the Cerro Tololo and European 
Southern observatories. 



3 MAKING MEASUREMENTS AND SUPPLEMENTARY 
STUDIES 

The observations with the MASS were made from August 2005 to November 2007. The device 
was attached to the photographic refractor telescope AFR-2 of 23 cm aperture and f/no 10. 
Valid data are only obtained after a proper adjustment of the Fabry lens which sets the sizes of 



effective wo rk apertures projected onto the telescope entrance aperture plane (jKornilov et al 



20031 . l2007h . The respective magnification factor equal to 15.4 was measured after installation 



and was confirmed by repetitive measurements in 2006. The physical sizes of the MASS segmen- 
tator and the effective apertures of respective channels are given in Table [TJ One can see that 
aperture size for the channel A is smaller than a typical Fresnel radius yXH = 3-7 cm depend- 
ing on turbulence height while the D-channel aperture is larger which is a must for successful 
application of profiles restoration algorithms. 

The weighting functions converting the turbulence distribution along the line of sight into 
scintillation indices depend on spectrum of light. Correct account of this factor is implemented 
using typical energy distributions in spectra of program stars and carefully the measured re- 
sponse curve of the device. While the device itself is well studied in a laboratory, the refractor 
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Figure 1: Colour equation of MASS at the refractor telescope AFR-2. Asterisks - real stars, a 
solid curve - parabolic approximation, circles - computed colour dependencies for a set of used 
typical spectral energy distributions after the applied response curve correction. 



Table 2: The device parameters p (non-poissonity) and r (non-linearity) by channels measured 
in 2006. 



Channel p a v r, ns o>, ns 



A 1.010 0.001 22.8 3.5 

B 1.012 0.001 21.4 1.4 

C 1.008 0.001 21.0 0.3 

D 1.007 0.001 19.4 0.2 



lens introduces significant absorption in the blue and hence strongly violates the actual response. 
Correction of this effect was made using photometric measurements of some 50 stars taken 15 
and 16 October 2006. 

Using the obtained atmospheric extinction coefficient (0 m .27 in the MASS bandpass), the 
colour equation mass — VvsB — V was constructed (Fig. [1]). The slope of this relation allowed 
us to correct the response of the system given an assumed typical glass selective absorption lav|l 
The original MASS response curve and those after correction are shown in Fig. [2j 

Unbiased estimation of turbulence strengths, especially when it is weak, is only possible 
given the correctly measured detector non-poissonity factors p in A and B channels. The flux 
in C and D apertures is normally strong enough to neglect this effect. Non-poissonity mea- 



surements involve the contro l light source which is built into the device (jKornilov et al 



v mea- 

1H; 



Tokovinin Sz Kornilov . 20071 ). Control measurements give also the non-linearity (dead-time) 



measure r which, contrary to p, is most important for C and D channels where bright stars may 
develop more than 10 6 counts per second. The results of these supplementary measurements 
made in two separate nights of 2006 are given in Table EJ 



x See Kornilov V., The verification of the MASS spectral response. 
http: //dragon, sai .msu. ru/mass/download/doc/mass_spectral_band_eng.pdf , 2006. 
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Figure 2: Original (black dots) and corrected (empty circles) spectral response curves of a 
MASS + telescope system. 



The statistics of observation time distribution by seasons is given in Table EJ A total of 
280 nights was occupied by measurements of which a small fraction (25 nights) were instances 
of short estimates of less than half an hour duration. Total measurement time is 1022 hours. 
The gap in June - September 2006 is due to device temporal failure. 

Inspite of significant gaps in measurement periods, the observation time distribution fol- 
lows more or less a typical seasonal clear nights allocation at Mt. Maidanak except for December 
and January when, by various reasons, no scintillation was measured. The pronounced grouping 
of measurements in the five periods is well seen: fall 2005, spring and autumn 2006 and the ones 
of 2007. 

An example of MASS measurement results is shown in Fig. [3] for a typical night of Septem- 
ber 9, 2005. The sporadic, usually quite short-term, increase of turbulence strength is seen in 
different layers (altitudes). The turbulence evolved from the evening time gradually declining 
in altitude. The end of night is dominated by the turbulence at 4 km. 

During the campaign the Maidanak observatory differential image motion monitor 
(DIMM) was mainly used in other sites of Uzbekistan. Simultaneous measurements with both 
devices were made in 2005 from 17 of September to 5 of November (18 nights), in 2006 - from 
11 of February to 31 of May (24 nights) and continuously during 24 nights from 25 of August 
to 16 of September of 2007. The results of MASS and DIMM data inter-comparison will be 
presented in a separate paper. 



4 BASIC MEASUREMENT RESULTS 

The MASS device is working under the Turbina program control ( Kornilov et al. . 20031 ) which 
operates under Linux OS. Although the functionality of this program includes not only the 
measurement control but also the real-time data processing spanning from scintillation indices 
calculation to the restoration of a current altitude turbulence profile, the repetitive off-line data 
reprocessing allows one to obtain more reliable and homogeneous results. This improvement is 
related to temporal drifts of the device parameters over 2-3 year period and, more importantly, 
to the continuous software update throughout the campaign time. 
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Figure 3: The September 9, 2005 results. Top: six fixed altitude layer representation of 
turbulence (see section H]). The stripe width is proportional to the layer power which scale is 
shown by the vertical bar to the right of the graph. Bottom: free atmosphere seeing j3f ree 
behavior (circles) during the night; crosses depict seeing formed above 2 km. 
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Table 3: Measurement duration in nights and hours by months of 2005 - 2007 



Season 
Months 


2005 
nights hours 


2006 
night hours 


2007 
night hours 


January 










1 


0.1 


February 






8 


20.5 


4 


15.5 


March 






12 


27.2 


6 


8.8 


April 






17 


71.6 


6 


15.7 


May 






23 


84.1 


10 


23.5 


June 










2 


4.8 


July 










26 


72.5 


August 


14 


58.8 






25 


87.8 


September 


30 


210.0 






21 


105.1 


October 


23 


73.8 


14 


15.9 


20 


61.4 


November 


5 


34.2 


5 


9.6 


7 


18.5 


December 















Such a reprocessing^] was performed with help of a special version of an Atmos program 
which is part of Turbina software, and a set of s/ie//-scripts which filter out the input data glitches 
and invalid output. Although observations were performed with data accumulation time for a 
single profile restoration equal to 20 s due to imperfect tracking of AFR-2, the data reprocessing 
has regrouped the input into commonly accepted 1 min accumulation time. 

As an outcome, more than 50 thousand integral atmosphere parameters sets and altitude 
profiles were computed. Each profile is quantified by turbulence intensities of six layers which 
are centered at 0.5, 1, 2, 4, 8 and 16 km. It should be stressed that provided values are not 
C\(h,\) but the integrated ones: Jj = C^(/ij)A/ij. Hence the input of each layer into the image 
quality is defined by the following relation: Pi = 2.0 x 10 7 • J^ 5 for wavelength A = 500 nm. 
The seeing of l'.'O corresponds to the integral turbulence strength J = 6.8 x 10~ 13 m 1//3 , and 
(3 = C'5 is produced by J = 2.14 x lO^m 1 / 3 . 

In Fig. [5] we show the cumulative distributions of the inputs of different atmosphere parts 
into the image quality. In other words, it shows the seeing for an observer at a given altitude 
above the site. Additionally the free atmosphere seeing (5f ree distribution is shown, which was 
calculated as integral atmosphere parameter without profile restoration. One can see that these 
distributions are statistically indistinguishable. 

The median layers contributions are given in Table H As a good comparison, the last 
column shows the median turbu lence distribution for a Cerro Tololo observatory quoted from a 
paper by Tokovinin et al ( 2003al ). 



An overall (3f ree measurement media n is / /46. A somewhat smaller value of 0'.'40 was 



obtained for the Cerro Pachon observatory (jTokovinin &: Travouillonl. l2006h in 2003 and 2005 



while those value for Las Campanas site carnp^n is 0"48 flThomas-Osin et all lioosh. 

Interestingly, the formal median value precision for such a voluminous data set is only 
0.01. Such an estimate of uncertainty is not adequate due to non-stationary nature of the 
considered phenomenon. The characteristic variation of a median seeing is illustrated by five 
seasonal medians. 



See Kornilov V., Shatsky N., MASS Data Reprocessing, 
http : //dragon. sai .msu.ru/mass/download/doc/remass .pdf , 2005 . 
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Figure 4: Cumulative distribution of the atmosphere layers input into the image quality: 1 - 
16 km layer, 2 — 8 km layer and above, 3 — 4 km and above, 4 — 2 km and above, 5 — 1 km 
and above, 6 — all atmosphere from 0.5 km and above (/3/ ree ). Crosses — (3 free distribution 
obtained directly from scintillation indices without profile restoration. 



5 ISOPLANATIC ANGLE AND ATMOSPHERIC TIME 
CONSTANT 



Isoplanatic an gle 9p is involve d to characterize a correlation of wavefront distortions in different 
directions (see Roddier . 1999). It is computed from scintillation data as an integral parameter. 
Cumulative 9q distribution is presented in Fig. [5] where the overall data distribution is accompa- 
nied by seasonal ones. The median overall measurement is 2'.'19 while it varies from 2'/35 in fall 
to l'.'92 in spring. Maximal median value was observe d in autumn 2007 — 2f.'4:7. As a reference 
it is worth to note the GSM es timat e of 1998 - 1999 dEhgamberdiev et all l2000l ) equal to 2'/48 
while Kornilov and Tokovinin ( 200ll ) measured isoplanatic angle within 2 / .'12 to 2 / /40 range. 

In practice, the most interesting is a 0q value in weak turbulence conditions. In Fig. [5] 



Table 4: Median contribution of atmosphere layers into the image quality (arcseconds) . 



Layer 


Median 


2005A 


2006S 


2006A 


2007S 


2007A 


CTIO'03 


0.5 km and above 


0.46 


0.43 


0.57 


0.41 


0.49 


0.43 


0.55 


1 km and above 


0.37 


0.35 


0.46 


0.32 


0.41 


0.34 




2 km and above 


0.32 


0.30 


0.41 


0.27 


0.35 


0.28 


0.43 


4 km and above 


0.29 


0.27 


0.38 


0.25 


0.32 


0.25 


0.37 


8 km and above 


0.20 


0.19 


0.28 


0.18 


0.22 


0.16 


0.29 


16 km and above 


0.14 


0.14 


0.15 


0.12 


0.14 


0.12 


0.16 
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the Oq distribution is shown for that half of measurements when seeing is less than median: 
Pfree < C/46. The median of such a scope of 9q is 2'.'57. 

It is natural that 8q is vastly determined by the upper atmosphere turbulence. In this 
respe ct its value should be typ ical for majority of sites located at 2.5 km elevation above the sea 
level (|Tokovinin et all 12003 j ). 




Figure 5: Cumulative isoplanatic angle 9q distribution for various seasons and conditions: solid 
line — total scope distribution, dashed - in springs, dot-dashed — falls, crosses — fall 2007, 
circles — isoplanatic angle when /3/ ree is better than median. 

Equally important for adaptive optics performance, apart from isoplanatic angle, is an 
atmospheric time constant (wi dely also called a coherence time) tq which is another by-product 
of scintillation measurement ( Tokovinin . 2002 ). The dedicated study show^l that the time 



constant from MASS estimation is systematically (factor of ~ 1.3) lower than observed by other 
methods. Meanwhile we present uncorrected values hereafter. 

Median of tq distribution shown in FiglB] equals to 3.94 ms for the total campaign and 
2.82 ms and 4.44 ms for springs and autumns, respectively. When restricted to Pf ree < C/46 
conditions one obtains a 5.41 ms median. If we apply an empiric correction mentioned above, we 
come to the time constant about 7 ms in good conditions. The latter coincides w ith the median 



time constant at Antarctic plateau Dome C obtained with MASS device in 2004 (jKenyon et al 
20061 ). 



One can recall the tq = 12 ms estimate of Ehgamberdiev et al (j2000l ) derived during the 
international astroclimate campaign of 1998 at Maidanak. Meanwhile, this value was based on 
a relatively short-term scope of measurements. 



3 See Tokovinin A., Calibration of MASS time- constant measurements, 
http: //www. ctio .noao . edu/^atokovin/prof iler/timeconst .pdf , 2006. 
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Figure 6: Cumulative distributions of atmospheric coherence time tq: spring seasons — dot- 
ted line, autumns — dash-dotted, full data scope distribution — solid. Black circles depict 
distribution when (3f ree is better than median. 

6 PARTICULAR TURBULENCE PROPERTIES OF FREE 
ATMOSPHERE ABOVE MT. MAIDANAK 



In order to study the role of different atmosphere layers, we build the dependence of the relative 
layers intensity Ji ayer / ' J total on the total turbulence strength Jtotai according to MASS measure- 
ments (i.e. excluding the ground layer). These ratios are shown in Fig. [7] being averaged over 
400 profiles. The layer input behavior differs strongly for the boundary layer (0.5 - 1 km) and 
the upper atmosphere (8 - 16 km). The former shows its input nearly proportional to the total 
power while the role of the upper atmosphere diminishes with the total growth of turbulence. 
Most representative here is the 16 km layer for which the relation is practically inversely propor- 
tional. This means that its intensity is almost invariable and absolutely insensitive t q proc esses 
in the boundary layer. Similar ph enome non is noticed by Tokovinin and Travouillon (l2006h and 
in a paper by Thomas-Osip et al ( 20081 ). 

On the opposite, the boundary turbulence is responsible for image degradation at Maid- 
anak. Keeping in mind the values J = 0.2 x 10 _12 m 1 / 3 equivalent to C//5 seeing, we are mostly 
interested in the graph beginnings up to values 0.5 x 10 _12 m 1 / 3 . At the 2 km height the turbu- 
lence intensity is seemingly lower than at other altitudes. 

As a support for this fact of a key role of the boundary layer, we present the Fig. [8] 
where the cumulative distributions of layers intensities are shown for a good seeing data subset 
(image quality better than median). The prominent feature is a nearly absent lower (1 and 2 km) 
turbulence for 50% of cases. The median of intensity in 4- and 8 km layers is 1.3-J-1.5 x 10~ 14 m 1//3 
which is slightly more than the median value for the lowest 0.5 km layer (1.1 x lO^^m 1 / 3 ). It is 
natural to think that at 0.5 km we see the development of the turbulence generated by a ground 
relief. 

The special case here is again the upper atmosphere at 16 km. The steep rise of the 
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Figure 7: Relative input of atmosphere layers in a total turbulence as a function of latter 



cumulative distribution tells about nearly constant turbulence which is quite significant. In half 
a time its intensity is in a narrow span of 1.5 -j- 2.8 x 10 _14 m 1 / 3 . Note that the high altitude 
turbulence is restored with a minimal error from the scintillation data. 



7 CONCLUSIONS 



An estimation of a total turbulence power from DIMM measurements ( Ilyasov et all 19991 ) allows 
us to evaluate the dominating input of the ground turbulence (up to 500 m above surface) to be 
65% on average and the first 200 - 300 m above surface give already about 50 % con tribution. 
The same contribution is found at Cerro Pachon by Tokovinin and Travouillon ( 20061 ). 

In cases of exceptionally good images (25% quaritile) the full image quality from DIMM 



measurements in 1996 to 1999 (jEhgamberdiev et all |200Cj) constitutes 0'.'55. Comparing with 
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data in Fig. 0] we come to the ground layer contribution of 60% in these conditions while the 
whole boundary layer inputs more than 70% leaving less than one third of a power to the rest 
of atmosphere. 

In such a situation the development of an adaptive optics system working in the visible 
should be directed to correction of the low turbulence. It is hardly possible to achieve the 
Strehl ratio improvement factor more than 5 here, meanw hile the large corrected AO field of 
view (of the order of several arc-minutes, see Iptigautl . [2002) will be a reasonable compensation. 
It is worth to note here that the maximum gain for the median Maidanak conditions (Fried 
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Figure 8: Cumulative distribution of turbulence intensity in 0.5 km layer (solid line), at 1 km 
(short dashes), 2 km (long dashes), 4 km (dots), 8 km (short dash-dots) and 16 km (long 
dash-dots) for cases of good images (f3f ree < C'46). 



radius r$ = 0.15m a t wav elength 500 nm) at the 1.5 m telescope AZT-22 computed according 
to Roddier formula ( 19981 ) for optimal atmospheric distortions correction is equal to 32. 

A profiting advantage of Maidanak observatory is a large time constant value to ~ 7ms. 
This lowers the critical demand of a system response time by a factor of nearly 2.5 and thus makes 
about three times more star s available for wavef ront referencing compared to Cerro Tololo and 
Cerro Pachon observatories ( Kenyon et al. . 20061 ). Such a good to also favours the development 
of optical interferometry at this site. 

This research was supported by Russian Basic Research Foundation (grant 06-02-16902a). 
The staff of Maidanak observatory helped a lot in organization and conduction of measurements 
for which we express gratitude to all of them. 
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